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Abstract 

A discretized massless wave equation in two dimensions, on an appropriately 
chosen square lattice, exactly reproduces the solutions of the correspond- 
ing continuous equations. We show that the reason for this exact solution 
property is the discrete analog of conformal invariance present in the model, 
and find more general field theories on a two-dimensional lattice that exactly 
solve their continuous limit equations. These theories describe in general 
non-linearly coupled bosonic and fermionic fields and are similar to the Wess- 

Zumino-Witten model. 
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I. INTRODUCTION 



Since our spacetime may ultimately prove to be discrete, field theories on discrete space- 
times have developed an interest in and of themselves In particular, field theories on 
two-dimensional discrete spacetimes can be interpreted as string theories with a discrete 
string worldsheet and a continuous target space. Exact solutions for such theories can of- 
ten be obtained even after quantization (see, for example, [0). Yet another approach to 
discretizing spacetime and constructing a theory of gravity based on stochastic properties 
of random lattices is advocated by R. Sorkin [^J. Investigation of field theories on discrete 
spacetime might therefore contribute to our understanding of discrete quantum picture of 
the Universe. 

We shall focus on the connection between the discrete-spacetime and the continuous- 
spacetime versions of the same field theory. This work was motivated by the (somewhat 
"mysterious" ) fact noted in Ref. |4[] that a particular discrete version of the two-dimensional 
wave equation gives exact solutions of the underlying continuous equation. In other words, 
the exact solution of the free scalar field theory coincides, on the lattice points, with the 
solution of the lattice version of the same theory. The main intent of this paper is to 
understand the origin of the exact equivalence of continuous field theories and their discrete 
counterparts, as well as to find more general classes of (classical) field theories that are 
exactly solved by properly chosen discretizations. The origin of the exact solution property 
of discretized equations turns out to be the discrete conformal symmetry of the models 
||. The exact solution property of the massless bosonic model has a direct application in 
numerical simulations of strings (which was the original context of Ref. |4]]). A quantized 
version of the discrete conformal symmetry may prove useful in string theory. 

The paper is organized as follows. We first present the exact solution property of the 
discretized two-dimensional wave equation. Then, in Sec. |HT| , an analysis of the discrete 
equations leads us to the concept of discrete conformal invariance (DCI), and we give general 
conditions for a given discrete field equation to possess DCI. In Sec. [TV] we present a general 
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description of discrete models with DCI, and show that they all yield exact solutions for the 
corresponding (first and/or second-order) differential equations; these solutions are always 
algebraically factorized. Further examples of equations with DCI based on Lie groups and 
semi-groups are presented in Sec. [V|. Conclusions follow in Sec. [VT|. 

II. EXACT SOLUTION PROPERTY OF THE WAVE EQUATION 

The two-dimensional wave equation came out of numerical modeling of cosmic strings 
J|] . The evolution equations for the target space variables X M (r, a) are 

X» T - X%, = 0, (1) 

with gauge conditions 

(X^±X£) 2 = 0. (2) 

Here, the subscripts r and a denote differentiation with respect to world-sheet coordinates. 

The string world-sheet can be discretized to a rectangular lattice with coordinates (r, o) 
and steps At, Act. A common second-order discretization of ([!]) would be 
X^ (r + At, a) - 2X^ (r, a) + X» (r - At, a) (r, a + Aa) - 2X» (r, a) + X^ (r, a - Ao) 



At 2 Aa 2 

(3) 

However, a special choice of discretization steps, At = Aa = A, was adopted in 0], and 
it was noted that the resulting discrete equations are not only simpler, 

X" (t + A, a) + X" (r - A, a) - X^ (r, a + A) - X^ (r, a - A) = 0, (4) 

[X"(r + A, a) -X"{r,a± A)] 2 = 0, (5) 

but in addition exactly solve the continuous equation ([I]). Namely, one can easily check that 
the general solution of ([[]), 
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X» (r, a) = a^ + (T-a) + at (r + a) , (6) 

(where a+ and at are vector functions of one argument, which are determined from the 
appropriate boundary conditions) exactly satisfies the discretized equation (|J). The gauge 
conditions (||) are then reduced to constraints on at, 

[a£ (x + A) - a£ (x - A)] 2 = 0. (7) 

To try to understand why this happens, let us examine the discretization in more 
detail. Expanding (|4|) in powers of the discretization step A, we notice that the discrete 
equation (01), on solutions of OX = 0, is satisfied to all orders in A: 



X (r + A, a) + X (r - A, a) - X (r, a + A) - X (r, a - A) 

d 2n X d 2n X\ 2 A 2 " 
dr 2n ~ da 2n J (2n)! 



^ x a2 ~ /a 2 "X 9 2n X\ 2A 2n 



n=2 

since 

a 2n x 9 2n x 



□ n x = o. 



Q T 2n Q a 2n 

Approximation to all orders means exact solution in the following sense. One can show that 

X« (r, a) = < (r - a) + (r + <r) , (9) 

where a+ and at are now arbitrary vector functions of one discrete argument, is in fact the 
general solution of (f|). Given a boundary value problem for the continuous equation ([!]), we 
can restrict the boundary conditions to the discrete boundary points and find the discrete 
functions a±. Then, from the form of the solution fl9|) it is clear that X^ (r, a) will coincide, 
on lattice points, with the solution of the original boundary value problem. This is what we 
call the exact solution property of the discretization (|J) . 

Like in the continuous case, the lines r ± a = const are characteristics of Eq. ([|). It is 
easily checked that the gauge conditions (0) are compatible with (^), in the sense that if the 
discretized gauge conditions @ hold at some point (r, a) , then the evolution equation (f|) 
guarantees that they hold everywhere along the characteristic that intersects the point (r, a). 



However, the exact solutions of the equations (0]) — © do not necessarily satisfy fl5|) for any 
discretization step A, because the continuous constraint (Q) imposed on the derivative of 
is non-linear and does not integrate to a discrete constraint (|7|). In this sense, the discrete 
constraints (|5]) provide only an approximate solution of (0). 

We have shown that a particular discretization (Q) of the wave equation exactly solves the 
corresponding continuum equation, regardless of the magnitude of the lattice step A. The 
remaining part of this work is intended to explain and explore this unexpected phenomenon. 
The "mysterious" cancellation of all higher-order terms in @ is due to conformal invariance, 
as will be explored in the next two sections. 

III. DISCRETE CONFORMAL INVARIANCE OF THE WAVE EQUATION 

Since the exact solution property is independent of the discretization step A, we are 
motivated to explore some kind of scale invariance in the model. We first notice that Eq. 
(|) does not mix the values of the field on "even" and "odd" sublattices, each sublattice 
consisting of points with r + a even or odd. Therefore we shall further consider only one of 
these two sublattices, which is itself a square lattice in the lightcone coordinates £± = r±cr. 
In Fig. 1 we drew four adjacent squares, each square corresponding to a group of four 
adjacent points on the lightcone lattice with the field values X^ 1 related by @, for instance 

X» (a') + X» (&') - X" (c') - X" (c) = 0. 

By adding together the corresponding equations of motion for all four squares in Fig. 1, it 
is straightforward to show that the same equation actually holds for points a, 6, c, and d 
that lie at the vertices of a larger 2x2 square, i.e. 

X M (a) + X» (h) - X» (c) - X^ (d) = 0. (10) 

In a similar fashion we find that the discrete equation (^) remains valid if we replace the 
discretization step A by its multiple nA, which means removing all but n-th points from the 
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lattice, and effectively stretching the lattice scale by a factor of n. We could call this fact 
scale invariance of the discrete equations. In fact, an even more general kind of invariance 
holds: namely, Eq. (Q) retains its form when applied to the four vertices of any rectangle on 
the lightcone lattice (i.e. with sides parallel to the null directions; such a rectangle would 
have its sides parallel to the £± axes). One can prove this by noticing that, by adding Eq. 
(U) applied to points a — e — a' — d and a 1 — d — b' — c, one obtains 

X" (a) + X» (b') - X" (e) - X" (c) = 0, 

which is the same equation applied to points a — b' — e' — c which lie at vertices of a rectangle 
rather than a square. The same procedure shows that Eq. (|j) applies to points f — a' — b — c 
as well. After deriving (|J) for the "elementary" 2x1 rectangles, it is clear that, for any 
given rectangle, one only needs to add together the relation (f|) for all squares inside it to 
obtain the same relation for the vertices of the rectangle. 

The fact that Eq. (^) applies to vertices of any lattice rectangle means that if we remove 
all points lying on any number of lines r + a = const and r — a = const from the lattice, 
the discrete equations will still apply to the remaining points. In lightcone coordinates, 
such a transformation amounts to replacing the (discrete) coordinates £± by functions of 
themselves: 

e = f + (e), r = r(r). (n) 

Here, the functions (£ ± ) are arbitrary monotonically increasing discrete-valued functions 
of one discrete argument. The monotonicity of these functions is necessary to preserve causal 
relations on the spacetime. 

We immediately notice a similarity between ( |TT| ) and conformal transformations in a 
two-dimensional pseudo-Euclidean space. As is well known, the wave equation ([I]) allows 
arbitrary conformal transformations of the world-sheet coordinates, 

f+ = /+(e+), r = r(r). (12) 
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This is the most general form of the coordinate transformation that preserves the pseudo- 
Euclidean metric 

ds 2 = dr 2 - da 2 = d£ + d£_ (13) 

up to a conformal factor. Transformations (|TT|) are obviously the discrete analog of conformal 
transformations ([12]) on the lightcone lattice. 

We have, therefore, found that the discrete wave equation (^) is invariant under the 
discrete conformal transformations (|ll|). We shall refer to this as the discrete conformal 
invariance (DCI) of the discrete wave equation. The natural question is then whether there 
exist other equations with the property of DCI, and whether such equations also deliver 
exact solutions of their continuous limits. In the remaining sections, we will answer this 
question in the positive, establishing a relation between the exact solution property and 
DCI. 

Note that the constraint equations (|5|) do not preserve their form under discrete conformal 
transformations. For instance, if the constraints hold for some step size A, 

[X" (r + A, a) - X* (r, a + A)] 2 = 0, [X 1 * (r, a + A) - X» (r - A, a + 2A)f = 0, 

it does not follow in general that they also hold for step size 2 A: 

[X" (r + A, a) - X" (r - A, a + 2A)f + 0. 

Accordingly, as we have seen in the previous Section, the discrete constraint equations do 
not exactly solve their continuous counterparts. 

IV. A GENERAL FORM OF A DCI FIELD THEORY 

To try to generalize Eq. (|J), we write the generic discrete evolution equation on a square 
lightcone lattice as 

X(t + A,<t) = F[X(t,o- + A),X(t,ct- A),X(t- A,a)}, (14) 
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where F is an unknown function. As we have found in the previous Section, the property of 



DCI will be satisfied if is invariant under the "elementary" conformal transformations, 



that is, if the relation (II) is valid when applied to the vertices of all 2 x 1 lattice rectangles. 



This requirement can be written as two functional conditions on F: 

F (a, F (b, c,d),b) = F (a, c, d) , (15a) 
F (F (a, c, d) , b, c) = F (a, b, d) . (15b) 

Here, a, b, c, and d are arbitrary field values which may be scalar or vector (or even belong 
to a non-linear manifold of a Lie group, as in our examples below), and F is a similarly 
valued function. Before we try to solve these conditions for F, we would like to show that 
for any function F satisfying flTo]), the discrete evolution equation ( |HD exactly solves its 
continuous limit. 

The continuous limit of fll4j) is obtained by expanding it in powers of A, for instance 

X (r, a + A) = X + X a A + ^ A 2 + O (A 3 ) , (16) 

and with the assumption that F (X , X , X ) = X , which is natural if we suppose that a 
constant function X = Xq must be a solution of (|l4|) , we obtain the following equations 
corresponding to first and second powers of A: 

X T = F 1 X a - F 2 X a - F 3 X T , (17) 

X TT = {F\ + F 2 ) X a(T + F 3 X TT + (F n + F 22 — 2F 12 ) X a X a + (F 23 — F 31 ) [X a X T + X T X a ) + F 33 X T X T 

(18) 

Here, Fi and Fij are derivatives of F with respect to its three numbered arguments (which 
are in general vector arguments, but we suppressed the indices in the above equations). One 
can easily verify that for F (a, 6, c) = a + b — c, these equations give the usual wave equation 

(I)- 

The derivatives and Fa are constrained by Eqs. flT5]). For example, the first derivatives 
satisfy 
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F\ — F\Fi, F 2 — F 2 F 2 , F 3 — — F\F 2 , (19) 

where Fi are understood as linear operators on the tangent target space. It follows that Fj 
are projection operators (i.e. operators P that obey P 2 = P) with eigenvalues and 1 only. 

Now we shall show that the equations (p)) — (P"8|) are exactly solved by ( fUj ) if the DCI 
conditions flI5| ) hold. Denote the exact solution of the continuous equations ([T7D — ( TB ) by 



X e (r,a), and the solution of the discrete equation (THf) by X d (r,a). Here, r and a are 



discrete lattice coordinates, and we assume that X e = Xd on the lines r ± a = (Fig. 2) . 



Since Eqs. (|17D — (JISQ were obtained from fll4|) by expansion in A up to third-order terms, 



X e satisfies the discrete equation (0) up to O (A 3 ), i.e. 

X d -X e \ r=2A ^ =0 = C(A)A\ C(A) = C + O(A), (20) 

in the limit of small A. Now we shall scale up the lattice by an arbitrarily chosen factor of 
n. As we have seen in the previous Section, from DCI it follows that Eq. (^4|) applies also 
to the vertices of an n x n square on the lightcone lattice, and this is a consequence of using 
Eq. flI3| ) n 2 times, once for each "elementary" square. At each elementary square, we can 
replace the field values X^ by X e and introduce an error of C (A) A 3 . An error of SX in 
X entails an error of not more than 3SX in F (X, X, X), because the first derivatives of F 
are operators with eigenvalues of and 1 only. Therefore, replacing Xd by X e in the n x n 
square entails an error in X (r = 2n, a = 0) of not more than 3n 2 C (A) A 3 : 

^-^e| T=2nAi(T= o<3n 2 C(A)A 3 . (21) 

However, we can also apply (|20|) directly to the n x n square to obtain 

^-^| T=2nAi(T= o = C(nA)n 3 A 3 . (22) 

It means that 

C(nA) < -C(A), 
n 

which forces C (A) = since C (A) is a polynomial in A. Therefore, the approximation 
is in fact exact. 
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This argument shows that it was in fact only necessary to expand ([14]) up to second order 
in A, and all higher-order terms will lead to differential equations which are consequences 



of (|PTD — (|i8|) , just as we have seen in the case of the wave equation. It also shows that the 
differential equations corresponding to given DCI field equations flU]) are always second- 
order or lower. 

Using the evolution equation (0), one can write the exact solution of the continuous 
equations for boundary conditions on the lightcone r ± a = 0. To find X (r, a) at a point 
(r, a) within the future lightcone of the origin, we construct a lattice that has (r, a) as one of 



/ T+a T+a] 




( T—a T—a 


\ 2 ' 2 J 


- 


V 2 ' 2 



and explicitly write X (r, a) through the boundary values: 

x M = F (* (I+£, I+£) , X ((IZ£, - IZ£)) , A" ,0, 0)) . (23) 

Now we try to describe a general class of functions F satisfying flTSp. First, we find by 
combining Eqs. ( |ToD that 

F (a, F (d, b,c),c) = F (F (a, d, c) , b, c) . (24) 

If we denote 

a* c b = F(a,b,c) , (25) 
the above will look like an associative law for a binary operation * c : 

a * c (d * c b) = (a * c d) * c b. (26) 
The condition F (a, a, a) = a looks like a unity law 

a * a a = a, (27) 



although it doesn't follow from ( 15 ) that a * a b = b for all b. However, if we suppose that 



is actually an operation of group multiplication, with the usual unity law 

a * a b = b for all b, (28) 
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and the inverse operation b 1 , then it is shown in the Appendix that there exists a 
reparametrization r of the field values such that the operation * c is written as 

r (a * c b) = r (a) * r (c^ 1 ) * r (b) , (29) 

where by * we denote the group multiplication in the group of field values. Note that 
although Eq. (|29"D is not symmetric with respect to interchange of a and b, such interchange 
is equivalent to a reparametrization r (a) = a~ l (where a -1 is the group inverse of a). 

We arrive at a picture of a field equation of the type flUJ), derived from a group multi- 
plication in an arbitrary Lie group. We shall assume that any needed reparametrization is 
already effected, and that the evolution is directly given by ( |I3D with 

F{a,b,c) =a*c- l *b. (30) 

For example, if we consider a vector space as an Abelian group with addition of vectors as 
the operation *, we again obtain the formula 

F (a,b,c) = a- c + b (31) 

for the discrete wave equation. 

In view of Eq. (30) , the exact solution (|23| ) has an algebraically factorized form, as a 



product of functions of the lightcone coordinates: 

X(r,a) = a(r + a) *b(r-a) (32) 

(the constant c" 1 is absorbed by either of the functions). This is a characteristic feature of 
the DCI equations we are concerned with. 

The continuous limit of the discrete equations fll4|)— (|30f) can also be described in terms 
of an arbitrary Lie group G as a target space; it is a Wess-Zumino-Witten (WZW)-type 
model [0]. If X is a field with values in a (matrix) group G defined on a two-dimensional 
manifold M, the action functional of the WZW model is defined by 

L — ^—( Tr( X- l d a X ■ X- l d a X] d 2 M + — [ e abc Tr ( X^d a X ■ X-%X ■ X^dcX) d 3 B, 
2A Jm v ' 2Aix Jb v ' 

(33) 
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where B is an auxiliary 3-dimensional manifold whose boundary is M, and Tr is the matrix 
trace operation. For a specific choice k — ±4, the equations of motion in the lightcone 
coordinates become 

d T (x~ 1 d±X) = 0, (34) 

with the general solution 

X £_) = X + • X_ (£_) or X_ (£_) • X + , (35) 

with X± being arbitrary G-valued functions. It is immediately seen that the solution ( |35|) 
is the same as (|32|), if we (naturally) choose * to be the group multiplication in G. 

V. EXAMPLES OF FIELD THEORIES WITH DCI 

In this Section we present some examples to illustrate the constructions of Sec. [TV] and 
to give a physical interpretation of the field theories arising from them. 



A. Interacting scalar fields 



As was noted in the previous section, we can obtain the discrete wave equation by using 
the general formula (|30D on a vector space considered as an additive group of vectors. Since 
all commutative Lie groups are locally isomorphic to a vector space, we should take a non- 
commutative group to find a less trivial example. The simplest non-commutative Lie group 
has two parameters and can be realized by matrices of the form 



{a, a} 

The composition law of the group is 



1 



(36) 



e" a 
1 



b 

V° 'J 



e a+p a + e a b 
1 



(37) 
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or, written more compactly, 

{a, a} {p, b} = {a + /3, a + e a b} . (38) 

The inverse element is given by {a, a} 1 = {—a, — e~ a a}. 

A calculation shows that the continuous limit of the discrete evolution equation 



with the composition law fl38D is (in lightcone coordinates £± 

d 2 a 



0, (39a) 



d 2 a da da 

Note that the non-commutativity of the group leads to asymmetry with respect to the spatial 
reflection (i.e. interchange of £ + and £_), but at the same time spatial reflection is equivalent 
to reparametrization a — > —a, a — > — aexp(— a) which corresponds to taking the inverse 
matrix to (06|) . 

The exact solution of (09|) is obtained from the general formula (|23|) : 



a (r, a) = a + (r + a) + a_ (r — cr) , (40) 
a (r, a) = a + (r + a) + e a+(r+<T) a_ (r - a) , (41) 

with arbitrary functions a± and a±. 

The equations (|39|) were obtained from a two-dimensional group and describe a pair of 
coupled scalar fields. More generally, one may start with an arbitrary n-dimensional non- 
commutative group G and construct the corresponding discrete and continuous equations 
describing n coupled scalar fields. The group structure will then be reflected in the cou- 
pling of the fields: for example, if the group G has a commutative subgroup H, then the 
corresponding parameters will satisfy a free wave equation. This can be easily seen from the 
example above: the elements of the form {a, a = 0} form a commutative subgroup, and the 
corresponding equation ( [39a| ) for a is a wave equation. 
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B. Fermionic fields: the discrete Dirac equation 



So far, we have been dealing with scalar fields. Now we shall consider the possibility 
of DCI equations describing fermions. The 2-dimensional massless Dirac equation for the 
two- component fermionic field is 



l a d a ^ = 0, 



(42) 





f° 


0^ 








7 = 

























where the corresponding Dirac matrices satisfy the usual relations of a Clifford algebra 
\l a -, 7 b } — ^Q ah an d can be chosen as 



(43) 



Here, the metric g ab in the lightcone coordinates £ a is g + = g + = 1, g ++ = g =0. 
With the choice (|43|) , Eq. (|^) becomes 

d + A = 0, d.ifj 2 = 0, (44) 

where tpx^ are left- and right-moving components of the field i)j (£ + , £ _ ). 
The corresponding lattice equations are 



(r + A, a) - ^ (r, a + A) =0, 
^ 2 (r + A, a) - V> 2 (r,<7- A) = 0. 



(45a) 
(45b) 



Their solution is 



^Ai = h {r + cr) , 1P2 = /2 (t - o") 



(46) 

(Here, /i^ are functions of 



which is also an exact solution of the continuous equations (f£ 
discrete argument determined by boundary conditions.) 

Since we again discovered a case of an exact solution, we naturally try to write the 
lattice equations ( £5 ) in the form (JTjJ) — (|5(]|) . This is possible if we define the multiplication 
operation * by 
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* 



















(47) 



(Note how the left- and right-handed components of the field propagate to the left and to 
the right of the multiplication sign.) Such a multiplication operation is associative, but does 
not allow a unity element and is not invertible, which would make it impossible to write 
a * c _1 * b as in (|30| ) . However, this operation has the property that a* x*b = a*b regardless 
of the value of x, and therefore we can disregard c _1 in (|30|). A matrix representation of 
this multiplication operation can be defined by 

I, , \ 



i V>i 





V 



1 o 

^2 



(48) 



Since the derivations of Sec. [TV] only use the associativity of the multiplication operation 
*, all our considerations apply also to cases where this operation does not have an inverse, 
such as in the case of semigroups ||. We shall be interested in a specific example of the semi- 
group structure represented in (pf ), where ipi and ^2 can be, in general, multi-component 
fields; we shall refer to such a structure as a "fermionic semigroup" . 



C. Coupled bosons and fermions 

Heuristically, a Lie group generates bosons and a fermionic semigroup generates fermions 
in DCI field theories. The direct product of a group and a semigroup would result in a theory 
describing uncoupled bosons and fermions. An example of a model containing coupled bosons 
and fermions can be obtained from a semigroup built as a semi-direct product of a group 
and a fermionic semigroup. A semi-direct (or "twisted") product of two semi-groups S and 
G can be defined if S acts on G, i.e. if for each s G 5 there is a map s : G — > G such that 

s (gi * 92) = s (gt) * s (g 2 ) (49) 
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and 



si(s 2 (g)) = (s 1 *s 2 )(g) 



(50) 



(Here, the multiplication denotes the respective semi-group operation in S or G, where 
appropriate.) The semi-direct product of S and G is the set of pairs {s,g} with the multi- 
plication defined by 



{si, gi} * {s 2 ,# 2 } = i s i * s 2,9i * «i (g 2 )} 



(51) 



One can easily check that this operation is associative. Of course, a group is also a semigroup, 
and the semi-direct product construction can be applied to two groups or to a group and a 
semigroup as well. 

The existence of an associative binary operation on the target space is really all we need 
to build a DCI field theory. We can obtain a generic theory of this kind containing both 
bosons and fermions by taking the fermionic semigroup S and some Lie group G on which 
S acts. Such pairs (S, G) can be constructed for arbitrary dimensions of S and G. We shall, 
to illustrate this construction, couple the two previous examples and arrive to a model with 
two interacting bosons and three fermions. 

To do this, we take S to be the simplest fermionic semigroup of fl48|). As G we choose a 
group like one represented by (p6|), but with two parameters ai 2' 



( a \ 

e ai a 2 



{a, ai, a 2 } 



1 
1 



(52) 



To couple bosons with fermions, we need to introduce some action of S on G. For example, 
we can multiply the column of parameters ai j2 of the group G by the right-moving part of 
the matrix s, 



(1 o\( ai \ ( ~ \ 



ifj 2 



v a 2/ 



a 1 



(53) 
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which defines the action of {■i/'i, ip 2 } on {a,ai,a 2 } as 

({a,ai,a 2 }) = {a,ai,i/} 2 ai} 



(54) 



One can check that the conditions (f49|) — (|50|) hold for such an action. The multiplication 
law of the resulting five-parametric semigroup is 

{ipi,ip 2 ; a,ai,a 2 } * {(j> u (f> 2 ; P, b x , b 2 } = {(j) 1 ,'4> 2 ;a + f3,at + e a b^a 2 + e a ^ 2 &i} (55) 

with a matrix realization 



{tpx, %l) 2 ; a, ax,a 2 } 



e a ai 
e a ip 2 a 2 
1 



V 



1 0! 





(56) 



The multiplication in this semigroup does not allow an inverse operation. Nevertheless, the 
problem with c _1 in the general formula ([30]) is circumvented because we have chosen the 
twisting of the semigroup and the group in such a way as to make the product a * c * b 
independent of the ipi t2 components of c. 

The continuous limit equations in this model are 

d + ^a = 0, 

= 0, d.ip 2 = 0, 
<9 + _ai = d + ciid-a, 
d + a 2 = ip 2 d + ai. 

As can be seen from the above, ^i j2 and a 2 are fermions, while a and a\ are bosons, 
a\ is coupled to a and a 2 is coupled to ip 2 and a\. (Here, we formally refer to the fields 
with first-order equations of motion as "fermionic" .) Similar models can be constructed for 
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a larger number of coupled bosonic and fermionic fields. Note that the field a 2 which was 
a boson in the {a,ai,a 2 } model, became a coupled fermion after we added the fermionic 
sector. 

VI. CONCLUSIONS 

We have explored the phenomenon of the exact solution of continuous equations by their 
discretizations. We formulated the property of discrete conformal invariance (DCI), and 
showed that any system of lattice equations possessing DCI delivers exact solutions of its 
continuous limit. In this sense, the conformal invariance is the cause of the exact solution 
property; the continuous limit equations must also be conformally invariant (although not 
all conformally invariant equations are exactly solved by any discretizations). We found a 
class of lattice equations, based on Lie group target space, with the property of DCI; their 
continuous limit corresponds to a theory of WZW bosons. We also found a more general 
class of theories based on semigroups describing scalar fields and fermions, which can be in 
general nonlinearly coupled to each other. 

In all these models, solutions to boundary value problems can be written explicitly 
(see Eq. (|23|)) using the multiplication law of the group or semigroup at hand. Expressed 
in this fashion through the boundary conditions on a lightcone, the solutions are always 
algebraically factorized. 

In case of the wave equation ([!]) with gauge constraints (H), it was found that the dis- 
cretized constraint equations @ are also exactly solved by certain solutions of the discretized 
wave equation, however they are not conformally invariant and, correspondingly, the exact 
solutions of the equations ([[]) — (H) do not necessarily satisfy (§) for any discretization step 
A. The "compatibility" of the discretized equations (f|) and (|) is perhaps due to the simple 
algebraic form of the solutions (|9]). 

The author is grateful to Alex Vilenkin for suggesting the problem and for comments on 
the manuscript, and to Itzhak Bars, Arvind Borde, Oleg Gleizer, Leonid Positselsky, and 
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Washington Taylor for helpful and inspiring discussions. 

APPENDIX 

Here we show that if the binary operation a * c b on target space V is not only associa- 
tive but is actually a group multiplication in some group G, then the target space can be 
reparametrized so that the operation * c becomes, in terms of the group multiplication *, 

a * c b = a * cT 1 * b. (57) 

By assumption, for each c G V there is a one-to-one map g c : V — > G from the target 
space to the group G such that 

9c (a * c b) = g c (a) * g c (b) . (58) 

The first condition of fll5|) can be written as 

a * b (b* d c) = a* d c, 

or, if we take of both parts, 

g b (a) * g b (b * d c) = g b g d l (g d (a) * g d (c)) . 

If we now denote g d (a) = x and g d (c) = y, where x and y are elements of G, this relation 
becomes 

dbg^ (x) * gtga 1 (g d (b) *y) = g b g d l (x*y). (59) 

This shows that ghg^ 1 is almost a group homomorphism, and we can make it one if we define 

h (b, d) (x) = g b g d x (g d (b) * x) , 

where h (b, d) is a map G —>■ G. After this ( |59"D becomes 

h (b, d) (x) * h (b, d) (y) = h (b, d) (x * y) 
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for all x,y G G. Now, obviously h (a, a) is the identity map, and h (a, b) h (b, c) = h (a, c) for 
all a, 5, c G V. This means that h (a, 6) can be expressed as 

h(a,b) = X(a) [X(b)]-\ 

where A (a) is an appropriately chosen homomorphism of G, and [A is the inverse of 

the homomorphism A (£>). For example, we could choose an arbitrary element a G V and 
define a map A : V — >HomG by 

A (a) = /i (a, ao) . 
Now, if we modify the function g b (a) by a A transformation: 

g b {a) = [A(6)] _1 g b (a) , 

then we obtain 

9c (b) * g b (a) = g c (a) , 
which similarly means that g b (a) is of the form 

g b (a) = r (6) * [r (a)] -1 , 

where r : V — > G is an appropriately chosen 1-to-l map, and [r (a)]" 1 is the group inverse 
of r (a). 

Finally, we can put the pieces together and find that 

F (a, b, c) = g; 1 (g c (a) * g c (b)) = r' 1 (r (a) * [r (c)]" 1 * r (6)) , (60) 
which is the desired result. 
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FIGURES 
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FIG. 1. The discrete evolution equation applies to any four vertices of a lattice rectangle. 




FIG. 2. Comparison of the lxl and n x n lattice squares. 
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